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Abstract 

We study the baryon density dependence of the vector meson spectrum using the D4/D6 
system together with the compact D4 baryon vertex. We find that the vector meson mass 
decreases almost linearly in density at low density for small quark mass, but saturates to a 
finite non-zero value for large density. We also compute the density dependence of the 77/ 
mass and the 77/ velocity. We find that in medium, our model is consistent with the GMOR 
relation up to a few times the normal nuclear density. We compare our hQCD predictions 
with predictions made based on hidden local gauge theory that is constructed to model QCD. 



1 Introduction 



If light-quark hadron masses are dynamically generated by the spontaneous breaking of chiral 
symmetry as it is commonly accepted in QCD, then it should be possible to "unbreak" the 
symmetry by dialing such external conditions as temperature and/or baryonic density. This 
means that the mass should change as the quark condensate - that is the order parameter of 
chiral symmetry - changes when the the external background such as temperature and/or bary- 
onic density is dialled. This phenomenon that reflects the vacuum change in strong-interaction 
physics is one of the main lines of research in nuclear and particle physics. 

In the case of temperature, the great progress in lattice QCD calculations together with 
recent experiments at relativistic heavy-ion colliders has enabled one to start understanding 
what is happening at high temperature. However the situation is totally different for dense 
baryonic matter such as one expects to encounter in compact stars. Because of the famous sign 
problem, lattice methods cannot access cold baryonic matter at high density. While there are 
a large number of models anchored on phenomenological or effective field theory approaches 
available in the literature, there are no known reliable tools to address the density that exceeds 
the normal nuclear matter density no — 0.16 fm~ 3 . Thus the density regime involved in compact 
stars is theoretically a wide-open area. On the experimental side, the forthcoming accelerators 
(such as FAIR/GSI, J-PARC etc.) will probe higher densities, but it will be sometime before 
results will become available to help guide the theoretical efforts. 

It is the objective of this paper to approach the hadron-mass problem by a holographic 
dual method along the line of recent developments in using gravity/gauge duality in top down 
approach [1, 2]. The advantage of the holographic approach is that it can handle - albeit in 
certain limits - the strong-coupling aspects of gauge theory relevant for dense matter in terms of 
weak coupling bulk gravity solutions. The logical procedure would then be the following. First 
one checks that the hQCD model so constructed correctly postdicts the observables probed by 
nuclear physics experiments, which means up to the normal nuclear matter density uq. It turns 
out that this is not an easy task as we will see. Once the first step is taken, then one could 
proceed to make predictions beyond the normal density no in as yet un-probed regime. The 
dense system we shall holographically realize is the confining D4 with the compact baryonic D4 
brane geometry discussed in [3, 4]. 

The reason for resorting to the D4D6 model is because it can accommodate the quark mass 
geometrically and we can discuss the condensate {^ip), that we call "chiral condensate" (CC). 
Note, however, that since the quark in this model is not a Weyl fermion at zero quark mass limit, 
chirality is not defined. Thus it is perhaps more proper to call it "scalar density" of fermion 
rather than chiral condensate given that the symmetry involved is not the axial symmetry but 
rahter the rotational symmetry of the brane embedding. Nevertheless one may think of it in 
terms of chiral symmetry [1, 5], because, based on it, one can still obtain the Gell-Mann-Oakes- 
Renner relation, which is one of the signposts of chiral symmetry in hadron physics. 

The plan of this paper is as follows. In Section 2, we describe what has been predicted 
in the boundary (gauge) sector and where the theory stands vis-a-vis with confrontation with 
nature. This will be done anchored on an effective field theory of strong interactions that is 
considered to at least partially represent QCD. It will be referred to as "BR scaling" standing 
for a circle of related ideas addressing the issue of chiral symmetry in dense medium. The 
setup of brane embeddings figuring in our holographic approach to dense matter is given in 
Section 3. In Section 4, we apply the formalism to in-medium spectra of the pseudoscalar and 
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transverse vector mesons that we are concerned with. Numerical analysis of the mass vs. density 
relation is made in Section 5. In Sections 6 and 7, we calculate the pion velocity and dispersion 
relation. Comparison with the predictions made by the class of approaches in line with Brown- 
Rho scaling is given in Section 8 with concluding remarks. Several details left out in the main 
text are relegated to Appendix. 

2 In-Medium Scaling In "QCD-Motivated" Approach 

In order to streamline our line of reasoning and specify what we are going to compare with, we 
first describe the presently available approach to the in-medium mass problem from a "QCD- 
motivated" angle 1 . As mentioned above, there are no model-independent gauge theory (QCD) 
tools for dense matter (except at the extreme - perhaps physically irrelevant - density regime 
where perturbative QCD is applicable). What we will compare with is a scaling relation that 
combines different lines of reasoning anchored on the notion of hidden local symmetry (HLS for 
short) [6] combined with chiral perturbation. Specifically the approach we will focus on is what 
we shall refer to as "Brown- Rho scaling approach" (BRSA for short) [7] . 

Originally, the BRSA to how hadron masses behave as the vacuum structure changes 
at increasing density was derived from the skyrmion model implemented with a scalar meson, 
dilaton, associated with the QCD trace anomaly [7]. It can be given a simpler interpretation in 
terms of the constituent quark model to which the skyrmion description should be equivalent 
in the large N c limit in matter-free space and presumably in medium. In the modern version 
based on large N c considerations [8], the constituent quark model should work not only for 
very low energy but also for intermediate energy in which deconfinement and spontaneously 
broken chiral symmetry co-exist. In the large iV c limit, the constituent quark mass wiq and 
the constituent axial-vector coupling constant are independent of N c , and baryons (mesons) 
are loosely bound three-quark (two-quark) states. In [8], the A c -independent mg is a given 
quantity in the matter-free vacuum, and there is nothing one can say from the analysis how 
niQ will behave in dense medium. The assumption we make here is that being dynamically 
generated by chiral symmetry breaking, rriQ will vanish when chiral symmetry is "unbroken," 
i.e., when the quark condensation (qq) — > 0. Suppose further that the condensate can be driven 
to drop as density is increased, going to zero (in the chiral limit) at a critical density n c . Thus 
rriQ will continue to drop as density increases. 

Now the constituent quark is a quasiparticle in the sense of many-body theory. If one 
applies the quasiparicle picture to AN C quark systems where A > 1 for nuclei and nuclear 
matter, the issue of how hadron masses behave in A > 1 systems would crucially depend on how 
the quasiparticle mass behaves in medium. Thus in the naivest scheme where the quasiparticle 
picture is assumed to apply not just to a single hadron but also to many-nucleon systems, the 
ratio of the masses in medium (with density n) denoted with the asterisk should scale 

m* M (n)/m M « m* B (n)/m B = $(n) (1) 

with the number of constituent quarks involved in the hadrons dropping out. Here the subscript 
M denotes light-quark mesons other than the pion and B light-quark baryons, i.e., nucleons. 

1 This is a bit misleading terminology if not a misnomer in that the circle of reasoning involved is not directly 
backed by QCD-proper calculations. It is based on various yet-to-be confirmed results inferred from effective field 
theories that are supposed to model QCD. We shall adhere to that terminology to distinguish it from hQCD 
approaches. 
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This is identical to the scaling proposed by Brown and Rho a long time ago [7] based on a 
dilaton-implemented skyrmion model. Even if the constituent quark mass is independent of iV c 
in the vacuum, it is not known how it behaves in N c in medium and how it depends on the quark 
condensate (qq) as the condensate goes toward zero in the chiral limit as the chiral restoration 
is approached. Typically 3>(n) = 1 — Cn/n with C = 0.1 ~ 0.3 [9]. One can offer a plausible 
argument that the scaling factor 3>(n) at low density is related to the pion decay constant in 
medium 

Hn) BR « (2) 

However it is not really known what happens at a density greater than that of normal nuclear 
matter, except near the critical density for chiral phase transition for which hidden local sym- 
metry theory [6] makes a simple prediction. In this (HLS) model - which provides a field theory 
support to BR scaling, local gauge symmetry plays an important role in giving the scaling 

m*(ri)/m p ~ (qq)*(n)/(qq) as n ->■ n c . (3) 

This would imply that the vector meson mass goes to zero at the chiral transition. To distinguish 
this from the BRS prediction, we shall call this "VM" prediction. 2 One of the motivations of 
this work is to test the relations (1), (2) and (3) using holographic QCD. 

In comparing the hQCD prediction that we are making with the BR/VM prediction 
(l)-(3), we should clarify how the scaling behaviors could be assessed. In order to be able 
to say which scaling is correct and which is not, one should know how the scalings appear in 
physical observables. For instance, in the BRS approach, the scaling figures in nuclear physics as 
parameters of the effective Lagrangian used in nuclear forces and in correlation functions. Only 
when the quasiparticle picture is a very good approximation will the scaling appear in a simple 
transparent form. Otherwise it will be compounded with complex many-body dynamics and it 
will not be an easy task to single out the "basic" quantity that is addressed here. With respect 
to the present situation in confronting nature, one can say that the scaling (1) is consistent 
with - though not confirmed by - a variety of nuclear observables but there is nothing yet to 
support or refute the scaling (3). To address this scaling in hQCD is the purpose of this paper. 
Ultimately it will be up to experiments to establish the correct scaling relation and confirm or 
falsify the predictions of our model. 

3 Computational setup: D4 and D6 and the Baryon vertex 

In this section, we detail our strategy in addressing dense matter. We will focus on the low 
temperature, highly dense system, holographically realized as the confining D4 with the compact 
baryonic D4 brane geometry [3, 4]. The non-supersymmetric, confining D4 brane background 
with Euclidean signature, is given by 

ds 2 = (|) (V^dx»dx» + f(U)dxl) + (?) (^ + U^l) 

= 9s (^y\ F 4 = 2 -^e 4 , f(U) = 1 - , R* = n 9s NJl (4) 

2 VM stands for the vector manifestation fixed point in hidden local symmetry theory [6]. 



3 



This background is related to the geometry with Lorentzian signature by the double Wick 
rotation 

x 4 i — > t, t i — > x 4 , U i — > U KK . (5) 

o c/ 1/2 

The Kaluza-Klein mass is defined as the compactified radius of the X4 direction: Mkk = 

And the bulk parameters (Ukk,9s, R) are related to the gauge theory parameters (Mkk, 9ym, 

A) as 

A 2 A/ 2 

g ' = 2WW U KK = -XM KK ll A 3 = ^, A = s ?- M iV c . (6) 



By introducing a dimensionless coordinate £ as j^pj = \jr) d£ 2 , the bulk geometry is rewritten 

3/2 3/2 2 

^ 2 =(^) (di 2 + dx 2 + f(.U)dxl) + (|) (d£ 2 + ^ 2 ), (7) 



and U, £ are related 



(_u_ \ 3/2 1 f / n 3/2 , /M 3/2 \ , = A-(^o)- 3 \ 2 = ^ 
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Table 1: The brane profile : the background D4, the compact D4 and the probe D6 



3.1 The compact D4 baryon vertex 

In confining phase, the baryon vertex is a D4 brane wrapping on S 4 , see Table 1. For the 
computational convenience, we expand the internal S 4 as R 1 x S 3 . The background metric (7) 
takes the form 



[A 3/2 



ds 2 = (^J (dt 2 + fdxl + dx 2 ) + R 3/2 VU + d6 2 + sin 2 9dfl 2 3 



The induced metric on the compact D4 brane is 

ds 2 m= f^ 3/ V + ^3/2^7 



R) 



f/2 



1 + 



de 2 + sm 2 e<ml 



(9) 



(10) 
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and £' = d^/dO. The DBI action for a single D4 brane on which N c fundamental string is 
attached is 

+ 2ita'F) + ^4 j A G (4) 
= t 4 J dtde sm 3 6 -^4- /3 (£ 2 + £' 2 ) - F 2 + 3A t = J dtC D4 , (11) 



where 



r 4 = n/£lzg s 



^jfK N c Ukk 



2 2/3 2 8 / 3 (2vr/ 2 )' M4 (2vr) 4 /5 



r 27ra' , ~ 27ra' _ . . 

A t = -— A t F = —F te . (12) 

Since it reduces the number of equations to be solved, we introduce the Legendre transformed 
'Hamiltonian' 

Ud4 = F ? 4 — Cd4 
dF 

= nj dj9yJu>y 3 (e + Z»)y/D(0)* + S m s e, (13) 
where the dimensionless displacement, T)(0), is defined by 

Mm - m , B M£ ,14) 

8F v/(" 2 -/^ /3 )« 2 + 5' 2 )-^ 

Then the equation of motion for the gauge field is 

d e D{9) = -3sin 3 #, (15) 
from which by integrating, we get the solution D(0) 

D(0) = 2(2u - 1) + 3(cos 0-- cos 3 6), (16) 

where the integration constant v determines the number of fundamental strings. Note that vN c 
strings are attached at the south pole (9 = 0) and (1 — v)N c strings at the north pole {9 = it). 
Here we set v = and then impose smooth boundary conditions (£'(0) = and £(0) = £o) at 
the south pole, which means there is no string at the south pole or equivalently all fundamental 
strings are attached at the north pole. Note that the numerical solutions of baryon D4 are 
parameterized by the initial value, £o- 

We considered a single compact D4 brane with A^ c fundamental strings attached at the 
north pole. The tensions of D4 and fundamental strings are balanced to have a stable config- 
uration. Assume that all fundamental strings are dangled at the same point 3 , a cusp at the 



3 If we take a D4-D6 system with a confining background and consider the baryon vertex as compact D4, both 
D4 and D6 get deformed due to the interaction between them through the fundamental string (Fl). When the 
strings and D-branes are in contact, the length of the Fl tends to be zero while the D-branes are deformed to 
replace them for the minimum energy configuration. Therefore the strings connecting D4 and D6 have zero length 
and D4 and D6 are in contact at a point. 
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north pole of D4, and the location of the cusp is U c . The force from the D4 brane tension can 
be obtained by varying the Hamiltonian of the D4 brane with respect to U c while keeping other 
variables constant; 



&H 



dU, 



X/l). ( ' ^ , 5 • (17) 

■ ■ " fix other values V 1 ~ & / 

where Tp = ^u^k * s ^ e t ens i° n of the fundamental strings. To have a stable configuration, the 
system needs something to support the compact D4 because the tension of fundamental string 
is always greater or equal to the D4 brane tension. That something is probe D6 brane. 

3.2 Probe D6 brane 

Now, we put the probe D6 brane where the other endpoints of fundamental strings are attached. 
The string endpoints can be understood as point charges on D6 brane with A t being the gauge 
potential that couples to this point source. Then the bulk metric (7) can be written as 

ds 2 = i-J (dt 2 + dx 2 + f(U)dx 2 ) + - I - J (dp 2 + p 2 dn 2 + dy 2 + y 2 d<p 2 ) , (18) 

where the D6 brane world volume coordinates are (t, x, p,0 a ). We assume that the only function 
that depends on p is y(p), with the transverse direction <f> set to zero. The induced metric on 
D6 brane is 

ds 2 m = (|) (dt 2 + dx 2 ) + (J) 2 [(1 + y 2 )dp 2 + p 2 dQ 2 ] , (19) 

where 

1 A/ 2 1 A 2 

i-WOp, R 3 -- 2j ^, o^t^j^t- v KK = -XM KK il (20) 

The DBI action for Nf D6 brane is 

S D a = J dtC D6 = -N f ne J e-^y/det(g + 27ra'F) 

= -r 6 J dtdpp 2 uf^f(l + y*)-F^ (21) 

where 

r 6 = N f p 6 v 3 n 2g -^l M6 = ^-L^, p = ^^p^ (22) 

To see more clearly the dimensionality of the Lagrangian and other quantities, we introduce the 
dimensionless coordinates (£,U,p,y) rescaled by £o,Ukk, 

U = U KK U, £ = £ol, y = t y, p = tof>. (23) 

Hereafter we shall drop the tilde in all dimensionless variables unless ambiguous. We define the 
constant of motion Q and its dimensionless partner Q from the equation of motion for F : 

9C D6 p 2 ufF ^ 



9F ^W)-^ 
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This Q represents the number of point sources (number of fundamental strings) and is 
related to Q by 

q - _ ^M KK 

^~ 2ira>T 6 ~ 4V39^ 

To compare our model with the real world, we identify the number of fundamental strings Q 
with the number of quarks. Then the baryon density ns is expressed in terms of Q 

N B Q 32 1 / 3 3 ~ / M KK V (0^ 

with the normal nuclear density given by n = 0.16 fm~ 3 = 1.28 x 10~ 3 GeV 3 . We note that 
from [1], the validity of the solution Eq.(19) is guaranteed if the following condition is satisfied, 

4 « i « 1. <P) 



(WkkMkk) 1 " (f\ > 1, 



And this condition is rewritten as 

A 1 

1 N 

X«! ^TTI » A ( 28 ) 

where we have used Eq. (20). It means A should be larger than 1 but less than iV c . Therefore if 
we take N c =3, the 't Hooft coupling A cannot be greater than 2.5 for the validity of the solution. 
The "Hamiltonian can be obtained by Legendre transformation; 

— b — Tp! <->D6 

OF 

= r 6 J dpV(p)y / l + y 2 . (29) 

The equation of motion for the embedding y(p) is 

y y 2 y dlogV dlogV 



1 + y 2 (1 + y 2 ) 2 1 + y 2 dp dy 
and the Gauss-law constraint gives 



+ z y . 2 a ta ^ = ' ( 3 °) 



F = y (31) 



2 , 48/3 



We can solve Eq. (30) if the boundary condition is given. In the Q / case, there are funda- 
mental strings which connect the baryonic D4 brane with the probe D6 brane. In this case, the 



4 This is not the energy of the system but a kind of on-shell action. Instead of solving the equation of motion 
for At, we use Eq. (24), the constant of motion, which is the solution of At. 
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tension of D4 and D6 branes should be balanced so the D6 embedding y(0) has a finite value 
so as to have an angle at the tip. This angle is determined by the force balance condition : the 
force at the cusp of the probe D6 brane is obtained as 



FD6 



D6 



dU r 



fix others 



27ra' 



(32) 



The whole configuration (compact D4+F1+D6) is stationary if there is a condition to balance 
the forces between them. The number of baryonic D4 is Q/N c since the flux from Nq strings 
should be canceled by a single D4 brane; the ratio between the number of fundamental strings 
Q and Nc is the number of baryonic D4. The condition which makes the system stationary 
(p = 0, y c = £ c , the subscript c denotes the position of the cusp) is: 



F D6 



Q_ 



F D 4, 



which is the "force balance condition" (FBC), and is simplified as 

. c 

y c = — • 

Vc 



(33) 



(34) 



With this, the value of y c that satisfies the FBC is uniquely determined for given values of 




Figure 1: Left:The embedding solution of baryonic D4 and probe D6 with 
10 -4 (thick, blue), 6.7(thin,green), 16.7(dashed,black). 



0.1 and Q 



and m q . The embedding solutions of D4 and D6 are drawn in Fig. 1 for a fixed y v {oo) with 
varying Q. It is not clear why three embeddings have the same y^, but we draw the figure only 
for small p because of the lack of the space. As we can see, at the cusp, y(p=0), the slope of D6 
brane is different for each different Q even for the same y^. 

4 Pseudoscalar and Transverse Vector Meson 

As was considered in [1], the mesons are identified as fluctuating fields on flavor brane. Since 
we have a non-trivial background, the fluctuations living on the background with a non-trivial 
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vacuum are 



y{p,x?) = y v {p) + e 5y(p,af), <f>{p, x") = e 54>{p, x") 

A t (p,x») = A t (p) + e6A t (p,a?), A x (p,x") = e 5A x (p,x») 

Sy(p,xn = e~ ikx Y(p), 5<t>{p,x») = e~ ik ^(p) 

5Mp,x») = e- ikx a t (p), 5A x (p,x») = e~ ikx a x {p) 

where e is the expansion parameter which is very small, and £ 2 = p 2 + y 2 



U 



1 C 



U K 4V3 £ 



2/3 



, Co 



Uk 
4V3' 



(35) 



(36) 



The embedding metric of the probe brane is 



ds 2 



^) U^dxV + ^^[(l + yl)dp 2 + p 
+e 2 (d a Syd b Sy + y 2 v d a 54>d b 5$) dx a dx b + 2 e y v (d a 5y)dpdx a 



2 dn 2 



(37) 



We plug the field configuration Eq.(35) into the DBI action Eq. (21) and expand it to e 2 order 
to get the linearized fluctuation. The fluctuating Lagrangian is 



Cy = r 6 



+■ 



+ 



w 



R 3 ...2/3 



Co 2& 



1 



w f( F 2 + w f\ 2wfy v y v '(-4F 2 + 3w^(l + y' 2 )) 
v L-Y — - —YY 1 

250 goCl 

F 2 k 2 + {k 2 -w 2 )wf)) Y 2 

2F 2 + (1 + yl)) (F 2 - 3u,f (1 + £)) (1 + f v )(Ay 2 v - p 2 ) 



+ (F 2 - cf(l + yl)) (2F 2 - 9 W f (1 + yl)) - 2uf F 2 {1 + y 2 )y 2 }y 2 
for the scalar field, 



(38) 



T 6: 



p 2 ir 



2^L 

for the longitudinal gauge field, 



f (l + y> 2 ) , R3k 2 wf(l + y > 2 ) 



9o 



Co 



e 



(39) 



Fp 2 



„q 2/3 , 8/3 , 



2w 1 l\ v {2F 2 -wf{l + y' 2 )) 



Co Ct 



9o 



goCv 



Yai 



(40) 



for the cross terms, 



T6 



8/3 . 2 , « 3 fc 2 F 2 u,f + (^_A; 2 K(1 + ^) 



Co 



£3 
Si; 



(41) 
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for the Goldstone mode, and 



T6 



2^ 



4/3 a , R 3 k 2 F 2 + {w 2 -k 2 )wf{l + y' 2 ) 



w 



+ 



So 



w 



2/3 



£3 



for the transverse gauge field where 



9o = wf(l + y' 2 )-F 2 



Q 2 + p^w 8 / 3 



(42) 



(43) 



Note all variables are dimensionless; we just dropped tilde for all variables. The equations of 
motion for the Goldstone mode and the transverse vector 5 are 



( P -^.T 



p 2 



4/3 / 
w + a x 



p 2 vl 

P 2 



(i + y >2 H 



w 



2 R 3 k 2 R 3 



4 8/3 
P u + 



Co £o Q2 +p 4 8/3 



w 2 x R 3 



k 2 R 3 



A 

P OJ. 



+ 

8/3 



& & Q2 + p 4 8/3 



(44) 



For convenience, we will not compute the longitudinal vector and the real scalar that are coupled 
to each other. What we will do is to find the dimensionless eigenvalue 



Mi 



-w Q 



Co 



9wi 



and M 



,R 3 



9wi 



a x 



-W~ 



(45) 



which is the thermal mass in dense medium when k s , k x =0. We impose the Neumann boundary 
condition at p = and the Dirichlett condition at p = oo for Eq. (44): \&'(0) = ->■ ^(oo) = to 
solve the eigenvalue problem. One might wonder why we cannot impose the Dirichlett condition 
at p = 0. We discarded this possibility since it is unnatural to fix the y-position of the connection 
point, which can easily vibrate. The other possibility is discussed in Appendix. And from the 
meson spectra for the light quark systems, i.e., rf and p mesons [13], we get the Kaluza-Klein 
scale Mkk = 1.039 GeV and the asymptotic distance of D4-D6 yoo- And the 't Hooft coupling 
constant A can be determined from the relation m q = y v (oo) M g* x . 



5 Mass vs Density 

The eigenvalue of Eq. (44) is obtained for varying density. For very small values of n/no, the 
mass drops slowly and then in the medium density region the mass drops almost linearly in 
density. Such a transition from slow to linear dropping takes place abruptly as m q decreases. 
For large density, the mass saturates to a finite value as seen in Fig. 2. The linear dropping 
regime is a consequence of transition from concave up to concave down of the mass as a function 
of density. Assuming the physical situation corresponds more or less to the linear dropping 
regime, we predict the mass shift of the transverse vector meson to be given by 

m v (n) = my(0)$(n) (46) 

5 Here we mean by transverse vector the SO(2) transverse vector. In medium, the Lorentz symmetry, SO(l,3), 
is broken, so we cannot say anything in the Lorentz covariant manner, but we still have spatial rotation SO(2). 
The unbroken SO (2) symmetry has transverse vectors A x ,A y when z is the direction of the wave propagation. 
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where 

$(n) » 1-0.18- - — . (47) 
A no 

How does this "prediction" compare with Brown-Rho scaling in nuclear processes? A 
variety of indications coming from analyses based on the scaling described in Section 2, i.e. 
Eq.(2), are that the vector mesons p and oj scale in the vicinity of n = no as 

$(n)~ 1—. — r with « » 1/4. (48) 

v ; l + y(n/n ) 1 V ' 

This would be roughly reproduced by our model prediction Eq. (47) if A is about 6. Given 
that our model is valid for A < N c , we find it to be disingenuous to fix N c = 3. In fact, there 
are no physical quantities in our model that depend on iV c separately. Rather, they depend on 
the combination of iV c and g s . We may therefore pick the parameters within the regime where 
the theory is valid. Furthermore, the saturation in the dropping mass to a finite value at high 
density predicted in this model (see Fig. 2) is seemingly at odds with the VM prediction Eq. (3) 
which may be closer to nature although it has not been verified experimentally. 
For the Goldstone boson, we get a similar result: 

/ 6 n \ 

77v(n) «"y(0) 1 - 0.20--— , (49) 
\ A n / 

showing that the Goldstone boson mass drops slightly faster than that of the vector meson in 
the linear regime. This behavior is not in agreement with what is expected in QCD. In QCD, 
the properties of the pion with small mass are still protected by chiral symmetry so that the 
mass remains more or less the same in medium as in the vacuum. This aspect will be further 
commented on in Section 8. 

One might attribute the discrepancies we have here to the nature of the probe limit of 
our system. At high enough density, the number of the fundamental strings are high and their 
gravitational effect could become significant so the background geometry could be quite different. 
This means that we cannot ignore the back reaction of the baryonic D4. Stated differently, what 
is happening is that after passing a critical density, the system becomes unstable but there will 
be a stable geometry with full back-reaction of the baryonic D4 and the flavor D6. This should 
correspond to "beyond quenched approximation" in lattice gauge theory. As is shown in [17], 
the gravity back reaction of charge can restore the possibility of phase transition driven by high 
density. 

One can now see the geometric reason why the mass has to decrease; since the effective 
tention \p) = tqV(p) in Eq. (29) increases with increasing Q, it costs energy to bend the 
brane. Note that the fluctuation wave function describes the deformation of the brane. With 
our boundary condition ^'(0) = the wave function becomes smoother for large Q and the 
energy for this less-bending configuration should be lowered. Therefore the large Q, the smaller 
is the meson mass. This feature is captured in Fig. 7 in Appendix B. 
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Figure 2: Left :The vector meson p mass with y 1) (oo)=0.301. Right: The Goldstone boson rj 
mass. Our result shows that the Goldstone boson mass drops slightly faster than that of the 
vector meson. The normal nuclear density no was fixed by taking A = 6. 



6 The "Pion" (rf) Velocity 

The pion velocity 6 is defined as the ratio between the temporal and the spatial pion decay 
constants, /*//*. From Eq. (41), we define the pion decay constant as a wave-function renor- 
malization factor of (j), 

ei r f 2 r 

C 4 d,4> =2 dtd z xd t 0(j) + J dtcPxd^d^ 

VsZq J Kiy/ga V&a Jo q) 

The numerical results obtained from these formulae are given in Fig. 3. One sees that the spatial 
pion decay constant decreases while the temporal decay constant increases as density increases. 
As a result, the pion velocity decreases to almost 0.6 at the normal nuclear density. This is 
similar to the result obtained in [18] from the Sakai-Sugimoto model[2]. Our present result 
shows more rapid drop of the pion velocity since the spatial part of the pion decay constant does 
not increase at all. In QCD-motivated effective theories, one expects that the temporal pion 
decay constant decreases, not increases, in density and hence there is again a difference from 
the gauge theory sector. This aspect will also discussed in Section 8. 



7 Density-Dependent Dispersion Relation 

The dispersion relation is the relation between the frequency w and the spatial momentum 
k. And the frequency is numerically computed as the eigenvalue of Eq. (44) by varying the 
parameter k 7 . Fig. 4 shows the dispersion relation of the pseudo-scalar (left) and the group 
velocity (right). At zero density the dispersion relation is w 2 — k 2 = M 2 . What we want to see is 

6 Here we use "pion" as the pseudo-Goldstone mode, not the real pion meson. In our system, pion corresponds 
to the rj' meson. In QCD, the rf is pushed up in mass by the SU(1)a anomaly which is the 1/N C effect. In the 
large N c limit, it is degenerate with the pion. 

We compute —w in Eq. (44) by varying the parameter k because of the sign convention —M — —w + k . 
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Figure 3: The pion decay constants are given as a function of density. The temporal (spatial) 
decay constant increases (decreases) from 17 MeV to 21(15) MeV at the normal nuclear density. 
The pion velocity is presented as a function of density in the right panel. The pion velocity 
decreases when density is increased; at the normal nuclear density, n = uq, the pion velocity 



is about 0.85. 



how this relation is modified in dense medium. Given the dispersion relation, the group velocity 
is calculated by In our approach, the dispersion relation has two interesting features. First 
there is a sudden decrease of group velocity in the intermediate k region. The group velocity 
is linear in k when k is less than the rest mass(0.77 GeV for p meson) and v g is a constant 
when k is large compared to the mass. But when k is larger than 1, the group velocity suddenly 
decreases and saturates to a certain value, Fig 4. Secondly, the group velocity is not 1 in our 
dense medium (compare the zero density case, red line in Fig. 4). It slowly decreases when 
baryon density is increased. 




Figure 4: Left panel: The dispersion relation for the transverse vector with y v (oo) =0.303. 
The eigenvalue w is plotted. The solid lines correspond to Q = (red), .2(green), l(skyblue), 
6(pink) as a function of k. Middle panel: ^ as a function of k. Right panel: The group velocity 

v g = dw }^ |fc->oo is given and for each line. 



13 



8 Remarks and Conclusion 



We have solved the embedding function of baryonic D4 and probe D6 with FBC (force balance 
condition) with the embedding parameterized by y v (oo) and Q. The density dependence of the 
meson mass obtained therefrom is studied for the same y v (oo). We find no qualitative difference 
between the Goldstone mode and the transverse vector mode: their mass drops substantially, at 
the normal nuclear density but does not vanish at high density. Thus chiral symmetry does not 
get restored at any density. Both are at odds with what one expects in QCD: The pion mass, 
protected by chiral symmetry (though lightly broken), should not change much by density and 
chiral symmetry should get restored at least in the chiral limit at large density. 

As mentioned in Introduction, whether the vector-meson mass can be taken as an order 
parameter of chiral symmetry to be probed in physical observables depends on whether local 
fields are relevant degrees of freedom at the phase change driven by density. In the vicinity of 
nuclear matter density where the quasiparticle picture should be applicable, hadron masses, both 
nucleon and meson, do seem to drop proportionally to the pion decay constant (more precisely 
the time component /*). The drop is roughly 20%. Our model agrees with the BRS/VM scenario 
when we take 't Hooft coupling A =6. The value of A needed to fit the mass drop near no is 
around 6, which exceeds 2.5, the bound for the validity of the scheme for N c =3. The number 
of colors N c , however, figures nowhere in our analysis and hence need not be fixed to any value. 
Hence taking A=6 does not spoil the validity of our model. 

It should be stressed, however, that it is not obvious that the system could be described in 
terms of quasiparticles in the vicinity of the phase change. The system could be totally "melted." 
Up to day, searches for the evidence of vanishing vector meson mass in experiments have not 
yet come out positive. Therefore our hQCD result that the vector meson mass cannot vanish 
is not ruled out by nature. It remains to see what happens when 1/N C and 1/A corrections are 
made to the model. 

Let us consider the pion velocity for which there are results in effective field theories in the 
QCD sector. Effective field theory encompassing the current algebra relations and spectral de- 
composition and using linear density approximation predicts that the Gell-Mann-Oakes-Renner 
relation - which is satisfied in the matter-free vacuum by the D4-D6 model used here - is given 
in medium in terms of the time component of the pion decay constant /* as [14] 

{fllU)\m%/m«f = {qqY/m (51) 

which implies that as (qq)* drops at increasing density, the /* will also decrease. The same 
line of arguments indicate that the pion mass is largely protected by chiral symmetry so that it 
changes little in medium. What we found in this hQCD model seems not to be following the 
QCD-anchored intuition. 

It is intriguing and perhaps significant, however, that in this hQCD model, the left hand 
side (LHS) and the right hand side (RHS) of Eq. (51) are balanced even though (f^/f-n-) 2 
and {qq)*/(qq) separately behave differently, as one can see in Figs. 3 and 5. In our model, 
(qq)*/(qq) ~ 0.99, » 21/17 and m*/m„ « 0.8 at n=n . So in Eq. (51) LHS « 0.95 

and RHS ps 0.99. The GMOR relation appears to hold even at finite density - as it does in 
matter-free space - in our hQCD model. So if there were any defect in one, the same would be 
present in the other. 

What can be happening with the BRS/VM can be understood more clearly in hot matter. 
Suppose that the relevant low-mass degrees of freedom near chiral restoration are the pions only. 
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Figure 5: The graphical presentation of the GMOR relation is given as a function of density. 



The dashed line is for (/*// 7r ) 2 (m*/m 7r ) 2 and the solid line for the fermion condensate, 



At no with 't Hooft constant A=6, 



(V>y>)(rao) _ 



(ipip)(n) 



=0.99.. 



In this case, one can use chiral effective Lagrangian field theory to compute near the critical 
temperature the pion decay constant /* - which is found to be equal to the axial susceptibility 
ua - which is not zero, whereas the space component tends to zero, hence finding v w — > [10]. 
This result is drastically modified when the vector meson mass drops to zero (in the chiral limit) 
as in BRS/VM [15]. In this case, both /* and /* drop to zero, with the v w going to 1. What 
happens is that the the zero-mass vector meson cancels the pion contribution in /* . The situation 
is not as clear in density as in temperature. However simulating dense matter with skyrmions - 
which is essentially the same as the BRSA - shows that v w behaves similarly in density [16]. 

The glaring differences between the the two - hQCD and QCD-motivated - approaches 
raise the big question as to where the defects are and how to remedy them. While the QCD- 
motivated approaches are found to be consistent with observations up to the density probed by 
experiments, going beyond that density remains unverified and hence could very well be incorrect 
due to the lack of reliable theoretical tools to handle strong-coupling regimes. On the other 
hand, the hQCD-based models - seemingly inconsistent with available observations in hadron 
physics - can more readily access the strong-coupling regime within certain approximations. 
Connecting the two approaches will likely provide the means to unravel the structure of the 
poorly understood cold dense matter system. To improve the situation for hQCD we first 
have to consider the back-reaction of the gravity to the probe brane, that is, the gluon-fermion 
interaction effect. The situation is quite similar in the Sakai-Sugimoto model. It is possible 
that this increasing tendency in both the pion decay constant and the quark condensate for 
increasing density is due to the large N c nature. In that case, we need to consider the finite -/V c 
effect, which is beyond the gravity approximation. The result on the pion velocity v n presents 
also a conflict between the dual descriptions. The result found in the previous section is that at 
increasing density, the pion velocity decreases from 1 in the vacuum to near zero at high density. 
Model calculations based on the BRS/VM modeling QCD indicate the contrary, that is, the v n 
tends to go back to 1 near the chiral transition point. 
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Appendix 



A WZ term 



Taking from Chapter 19 of [12] the type-IIA super gravity solution of D4 brane 

1 



ds z 



VHd4 



[-Wdt 2 + dx 2 ] + v^D4 



(5) 

ty 1 y 2 y 3 y 4 



ae-*>(fl^-l) 



Hm = i + ^i, w = i + hD ^; a2 \ hm = 9 -^l, 
p 6 p 60J4 



(Al) 



the near horizon limit H~D4 = and renaming the coordinates, we write the non-extremal 
solution as 



ds 2 




[-fdt 2 + dx 2 




t i U k MS) 
where the parameters are identified as 



dlP 
[ f 
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+ u 2 dnj 



3\ 1/4 
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9s 
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R 3 
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p = U, Hda = R 
The five form field has its six form field strength, dC^ 



773 

9s, a 2 = l+ UK 



flDA 



dC {5) 

tx 1 x 2 x 3 x 4 U 



Ui u 2 



and its magnetic dual is 

0(41 - ^ c( 1^, = ^vyyw"<% w 



u 



1 R " sin 3 Q 3 ^1 + ^3^ = 3^ sin 3 6 n 3 A / 1 + 



C/ 3 



Note that the volume of S" 1 is 

d Sn -iV = sin n ~ 2 (0i)sin ri - 3 (02)---sin(0 n _ 2 )#i ■••#„_! 
d0 4 = sin 3 0cZQ 3 . 
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(A3) 



(A4) 



(A5) 



(A6) 



16 



B Other boundary condition 



In our treatment given in the main text, we used the Neumann boundary condition (hereafter 
N). Here let us consider the Dirichlet boundary condition (hereafter D) in computing the meson 
spectrum 

iJ/(0) = — ► *(oo) = 0. (Bl) 

The spectrum is the same as with the Neumann boundary condition at Q=0 but different 
at finite Q. The ground-state and the excitation spectra are presented in Fig. 6 with two 
different boundary conditions and y v (oo) =0.1 (left), l(right). The mass computed with D 
increases when the density is increased, while the mass with N decreases. The origin of this 
difference can be seen Fig. 7. With N, the wave function is more flattened for larger density, 
so the energy cost is lowered. But with D, the wave function bends increasingly more near 
the origin when Q becomes large. With D, the system costs more energy to fluctuate from the 
vacuum and with N, the energy cost is lowered as the density is increased. 



20 25 



Figure 6: The meson masses is presented as a function of the dimensionless density, Q. The red 
thick line is of the Dirichelet boundary condition and the blue dashed line is of the Neumann. 
TOP - The masses of lowest and first two excitations of r]f are plotted from bottom. Top-left 
figure is of yoo=0.1 and top-right is of i/oo = 1. BOTTOM - The masses of lowest and first two 
excitations of transverse vector meson are plotted from bottom and the same parameters are 
used as top figure. 
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Figure 7: Top : The wave function of <5(/>with Dirichlet (Left) boundary condition and with 
Neumann (Right) boundary condition. Bottom : The wave function of 5(f). 
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